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RESIN SELECTION CRITERIA FOR "TOUGH" COMPOSITE STRUCTURES 


by 

C. C. Chamis and G. T. Smith 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

Resin selection criteria are derived using a structured methodology con- 
sisting of an upward integrated mechanistic theory and its inverse (top-down 
structured theory). These criteria are expressed in a "criteria selection 
space" which can be used to identify resin bulk properties for improved com- 
posite "toughness". The resin selection criteria correlate with a variety of 
experimental data including laminate strength, elevated temperature effects 
and impact resistance. 


INTRODUCTION 

It is well known in the fiber composite community chat resins (matrices) 
provide the composite with the capability to resist load by xeeping the fibers 
in place. The capability of the resin to keep the fibers in place results 
from a combination of chemical, thermal, and mechanical interactions. These 
combined interactions produce the in-situ resin physical, hygral, thermal, and 
mechanical properties which provide the composite with the requisite struc- 
tural integrity in service environments, in general, and in turbine engine 
service environments in particular. A structured methodology is needed which 
can be used, in a formal way, to identify bulk (neat) resin characteristics 
which translate to quantifiable composite structural /mechanical behavior. The 
result of such structured methodology will be a set of criteria (guidelines) 
which can be used in advance to screen and/or select resins with the desirable 
bulk properties in order to provide the specified composite properties. The 
objective of the proposed paper is to describe a new structured methodology 
developed at Lewis for assessing, evaluating and identifying desirable bulk 
resin characteristics for specified structural composite integrity (fatigue 
resistance, fracture toughness, impact resistance, compressive strength, buck- 
ling resistance, vibration frequencies, and "toughness"). 

The structured methodology is based on an upward integrated mechanistic 
theory consisting of composite micromechanics, composite macromechanics, 
laminate theory and structural/stress analyses and its inverse (a top-down 
structured theory). All these are used formally to identify the resin charac- 
teristics which have a significant effect on composite mechanical behavior. 

The structured methodology is developed by using mainly matrix notation. Ex- 
panded equations are summarized in figures with appropriate schematics to 
illustrate the simulation and define the notation. The notation is also sum- 
marized in the Appendix for convenience. The structured methodology is based 
on Lewis' research activities on this subject during the past decade (for 
example references I to 3). The references cited in the text mainly refer to 
Lewis' research. However, these references include relevant references from 



the literature. The results obtained are summarized in convenient criteria as 
simple equation or ratio form which can be used to select resins a priori for 
improved and/or specific composite toughness. 


STRUCTURED METHODOLOGY 

The structured methodology used to develop the resin selection criteria 
embodies composite micromechanics, composite macromechanics, combined stress 
failure criteria, laminate theory and structural /stress analysis. The struc- 
tured methodology is integrated from composite mechanics upward to structural 
analysis in order to relate the structural response to constituent materials 
(fibers and matrix). It is then used in reverse, top-down, from structural to 
composite micromechanics. The structured methodology is described herein 
using a top-down approach in order to formally relate a specific structural 
response (for example displacement, stress intensity, stress wave propagation, 
impact resistance) to constituent material properties (fibers and resins). 
Resin and matrix are used interchangeably throughout this discussion. 


Structural Response 

The mathematical model describing the general structural response of a 
structure is given, using matrix notation, by 


[M]{ u> + CC]c 0} + [K]{ u> = { F ( t )} 


( 1 ) 


where {u>, (u> and { u> are the acceleration, velocity, and displacement 
vectors, respectively; [M], [C], and [K] are the mass, structural damping, 
and stiffness matrices, respectively; and { F ( t )} is the time dependent 
force vector. The natural frequencies and buckling resistance (buckling load) 
of a structure are described by special cases of equation (1), respectively. 


<tK] - »2[M]> 

= [oj 

(2) 

< [KJ - »2[i]> 

- [0] 

(3) 


where <o is the structure's natural frequency associated with a specific 
vibration mode, \ is tne eigenvalue containing the buckling load of a speci- 
fic buckled shape and [I] is the identity or unit matrix. If it is further 
assumed that [C] is proportional to [K], ([C] = y[K]), prior to any 
damage, then equation (1) can be written thus 


+ y[K]{ Cl} + [K]{ u> = { F(t )> 


( 4 ) 


Equations (2) and (4) depend on composite material properties embodied 
primarily in [M] and [Kj while equation (3) depends only on [KJ. For 
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most structural epoxies the density is about the same. 4 Therefore, the 
major influence of the resin in the global structural response of a composite 
structure is through the stiffness matrix [K]. 

The global structural stiffness matrix [K] is an assemblage of local 
stiffnesses. The local stiffnesses are related to the force deformation rela- 
tionships and the local geometry. The force deformation relationships for a 
composite laminate, including hygral (moisture) and thermal forces, are given 

Oyt 


{N> 

1 _ 

[AJ [C] 


<V 

| j tN T*| | 

<v 

{M} 


[c] T [ 0 ] 

- m 


{u} 

1 l< M T )i i 

‘V 


where { N} and {M} are the resultant force and moment vectors at the 
section, respectively. The subscripts T and M denote thermal and moisture 
forces. The reference plane strain vector is given by {e 0 } and the local 
curvatures by { u> . The arrays [A], [C], and [D] denote axial, coupled and 
bending stiffness, respectively. The various arrays on the right side of equation 
(5) for an N* ply laminate, using laminate theory^, are given by 


[A], [C], [D] = 


■ - V j <4 - z b> j - 2 b>] , 


{ Ny} , {My} 


IN 

■ E <([< z t -vi< z t 2 - z b>] *™ t ee] 

i = 1 XL 


( 6 ) 




n 

- L <Cf z t - W < z t - 

i = 1 


Am[R]'[E] { B} 




In equation (6), Z locates the i t(l ply through the thickness from a reference 
plane, [R] is the strain transformation matrix which defines the i th ply 
orientation (ply material axes) relative to the structural axes, [E] defines 
the ply stress-strain relationships, {a} the ply thermal expansion coeffi- 
cients, and {g} the ply moisture expansion coefficients. The relative 
changes aT and am denote the changes in temperature and moisture, respec- 
tively, of the i th ply measured from a reference condition. The temperature 
and the moisture for the i™ ply are determined by heat transfer and moisture 
diffusion analyses. 
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The important point to note from equation (5) is that the resin properties 
influencing structural response are reflected through the ply property arrays 
[E], {a}, and {6}. Another important poi^t to note is that the top-down 
structured theory described by equations (1) to (6) predicts global structural 
response which is in very good agreement with experimental data. For ex- 
ample: 1) the natural frequencies and mode shapes predicted by equation (2) 

are in excellent agreement with experimental data for fiber composite fan 
blades 6 and hybrid composite blades 7 ; 2) the buckling loads of anisotropic 
plates, including bending, stretching, and coupling, predicted by equation (3) 
are in good agreement with experimental data 6 *”; 3} the impact displacements 
predicted using equation (4) are in good agreement with the high-speed movie 
data for a large hybrid composite fan blade 7 ; and 4) the hygrothermomechani- 
cal response of a variety of angleplied laminates predicted by equation (5) 
are in very good agreement with measured data^. This good agreement of the 
predicted various structural responses with experimental data verifies that 
the global structural response of composite structures is formally related to 
the ply properties [E], fa), and {6} in equation (6). These ply properties 
can be formally related to matrix properties by continuing the top-down 
structured theory through composite macromechanics, combined-stress failure 
criteria, and composite micromechanics as will be described later. 


Stress Intensity/Concentration 

Experimental data from composite laminates with circular holes and slit 
type defects exhibit the same fracture characteristics. These characteristics 
are similar for tensile or compressive loads and for a variety of hygrothermal 
environments^- 6 . Typical results are shown in figure 1. Explicit equations 
are available which describe the stress concentration/intensity in the 
vicinity of a circular hole in an infinite composite angleplied laminate^. 

The equations, relevant to this top-down structured theory, are summarized in 
figure 2 with accompanying schematics. 

It can be seen in figure 2 that the stress concentration ratio ( CT c ©e /°cx X ^ 
in the vicinity of the hole depends on the laminate properties Ec ee * E C xx» ^cyy* 
Gcxy > and v c ^y. These properties are det<rmined using the following laminate J 
theory equation 6 


[E ] ■ I X) < t CR] T [EJLR j > 1 (7) 

c 1 . 1 


where t c is the laminate thickness and t-j is the ply thickness. Comparing 
equations (?) and (6), it can be seen by inspection that [E c ] = |/]/t c . 
Therefore, the stress concentrations in the laminate (composite-structure) are 
formally related to ply properties through equation (7) and the equations in 
figure 2. The hygrothermal effects on the stress concentration are determined 
from equation (5) since the local laminate stresses are related to {N> and 
{M}, or {e 0 > and { \i) as will be described later. 


Stress Wave Propagation 
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The speed (C) of the inplane stress wave propagation in an angleplied 
composite laminate in the X, Y, and XY (shear) directions, respectively, are 
given by 


C 

x 


9E 


cxx 


v cxy v cyx^ 


( 6 ) 


C 


y 



g£ 


SEL 


V V 

cyx cxy J 


1 

1 


(9) 




1 

2 


( 10 ) 


where g is the gravity acceleration, p is the laminate density, G is the 
shear modulus and the subscripts denote respective directions. The corres- 
ponding through-the-thickness normal and shear speeds, respectively, are given 
by 


z 

( gE C Z?J 

(ID 

\z 

(gG cx ,/ p ) 

(12) 

y* 

(gG cyz / P ) 

(13) 


The speeds predicted by equations (8) to (13), except equation (11), have 
been correlated with experimental data 1 . These equations have also been 
used extensively to theoretically investigate stress wave propagation in 
angleplied laminates due tc normal, oblique and edge impacts 1 **. 

Theoretical predictions of stress wave propagation due to point impact 
require a contact law in general. One such law for an impacting elastic 
sphere is given by the following approximate relationships 14 . 


F = (4/3)(R s ) 1/2 


E E 
s czz 


3/2 


E s * (l-v s )E c 


zz 


(14) 


where F is the contact force, R is the radius of the impacting sphere and 
6 is the local indentation. The subscript s denotes impacting sphere 
properties. 
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Taken collectively, equations (8) to (14) show that stress wave propaga- 
tion in composite laminates depends on the laminate properties E, G, v, and 
p. These laminate properties are formally related to ply properties through 
equation (6) as was already mentioned. 


Ply Strains and Stresses 

In the top-down structured theory, the i th ply strains { e-j > are formally 
related to global variables {e 0 } and i w> by the following matrix equation. 


{ ei ) - [Ri ] (c 0 > ~ Z i 


(15) 


where the glotal variables are determined from structural analysis. The cor- 
responding ply stresses are given by 


i o -i } = [Ei] i e) - am{g} - aT{q} 


(16) 


For the special case of in-plane loads only, the ply stresses are formally 
related to laminate stresses {o c > by 

{ oil = [Ei ][Ri ][E C ] 1 { CT c 1 (17) 


where [E c ] is given by equation (7). Equation (17) is significant because 
it shows that the ply stresses depend on local ply properties [Ei] as well 
as integrated laminate properties [E c ]. Equation (17) constitutes the ply 
stress influence coefficients for the case of in-plane loads only. 

The ply strains and stresses determined from equations (15) and (16), or 
(17), are normally comnared to ply fracture properties using available failure 
criteria such as ply: 


maximum strain: 

t ^i } < 

u t ) 

(18) 

maximum stress: 

i°i 1 1 

{5,} 

(19) 


or combined-failure stress criterion, for example ^ 5 


F(o.,S ,K ) = 1 
i i v 



„ a ilia? i22q 

tU \ 11 .^ 228 . 


( 20 ) 


where oi denotes the i th ply stress along the material axes denoted by 
the numerical subscripts and with sense denoted by a or 6 (tension or 
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compression). The corresponding ply uniaxial strengths (fracture stresses) 
are denoted by S. The coupling coefficient 12 depends on the ply 
elastic properties E all , E t 22. v t i 2 , and v t 23- 

The ply elastic constants [E £ ] and the uniaxial strengths (S t ) can 
be formally related to matrix material properties using composite micro- 
mechanics as will be described subsequently. 

Equations (18), (19), and/or (20) are used to assess laminate/composite 
structural integrity, durability, and/or composite toughness. This may be 
stated as the magnitude of stress {o c }, resulting from stress concentra- 
tion due to impact or defects, that a composite can sustain prior to ply or 
interpiy damage which will degrade (1) the composite global structural 
response or (2) the composite life/durability^. The composite life, or 
durability, is usually measured by its resistance to cyclic (fatigue), 
mechanical, hygral, and/or thermal loading. 


Composite Micromechanics 

Composite micromechanics is the discipline which formally relates ply 
properties to constituent material properties 4 . The properties pertinent to 
the development of "resin selection criteria" pursued herein are the ply 
mechanical properties (elastic [E A ] and strength S*) and the ply hygrothermal 
properties (b* and c^j where the subscript a (instead of i) has been used 
to denote ply properties, in general. In addition, the hygrothermal degrada- 
tion effects on the mechanical and thermal properties are related to the 
"hot-wet" ply environment and to the "hot-wet" glass transition temperature of 
the resin through a hygrothermal degradation factor Fm (HGTM). 

The equations for predicting ply properties in terms of constituents are 
summarized in figure 3. The notation in these equations corresponds to the 
schematic in the figure. The ply elastic properties are explicitly related to 
the matrix (resin) (r) and fiber (f) properties and to the hygrothermal 
degradation through (Fm). It can be readily observed in these equations that 
the "resin-controlled" ply properties are: E»22* ^i33» 6 a 12» ancl ^a23* 

The equations for ply longitudinal strengths (tension and compression) with 
attendant schematics for various fracture modes are summarized in figure 4. 

The dependence of the ply longitudinal compression strength on resin is 
through the resin properties (denoted by subscript (r)) and the ply intra- 
laminar shear strength (S i i2S)* The hygrothermal degradation effects on 
ply longitudinal strengths are incorporated through the respective resin prop- 
erties using the degradation factor (Fm) as shown in figure 3. Similarly, the 
equations for ply transverse (tension and compression) and intralaminar shear 
properties are summarized in figure 5. The equation for lower bound strength 
is derived assuming a resin slab (plate) perforated with a regular array of 
holes. The correlation between this equation and experimental data on yield 
stress of a steel plate^ is illustrated in figure 6, where S p and S 
denote yield stresses of the plate with and without the perforations, respec- 
tively. The parameter R p * denotes perforated area ratio which equals 
Kf or the fiber volume ratio. The corresponding equations for ply hygral 
{ 8^ ) and thermal (« t ) properties are summarized in figure 7. 

The equations summarized in figures 4, 5, and 6 close the loop in the 
structured methodology required to develop the "resin selection criteria". 
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Recall that this methodology consists of two multilevel theories: (1J An 

upward integrated theory which formally integrates resin elastic, hygral and 
thermal properties into composite structural response and (2) a top-down 
structured theory which formally relates the composite structural integrity/ 
durability to the same resin and, in addition, strength properties. The 
hygrothermal environment and degradation effects are included In both 
theories. 


RESIN SELECTION CRITERIA— DERIVATIONS AND IDENTIFICATIONS 

The resin selection criteria are derived with the aid of equations (17) 
and (20). Equation (17) is used in the form of Ply Stress Influence 
Coefficients (PSIC) defined as 


I >,] - [EJ[R,][EcH (21) 


where the subscript t is again used to denote ply properties in general. 

The expanded form of equation (2) when the composite is subjected only to 
°cxx 1S shown in figure 8. Analogous equations can be derived when the 
composite is subjected to other stresses or combinations*'. Using the PSIC 
from the equations (fig. 8), equation (20), accounting for hygrothermal degra- 
dation effects {f m ) and hygrothermomechanical degradation effects (/hjm). 
and rearranging results in the equation summarized in figure 9. The definitions 
for and ^TM are included at the bottom of the figure for complete- 
ness. The equation in this figure explicitly relates ply combined stress 
failure/damage (first ply failure) to the composite stress and includes the 
hygrothermomechanical degradation effects. The resin influence is through the 
properties with subscripts 22 and 12 and can now be readily assessed by 
expressing these properties in terms of constituent properties using the 
appropriate composite micromechanics equations from figures 3, 4, and 5. 

The results obtained for transverse strength, transverse modulus, energy 
density and first ply failure, respectively, are summarized in figure 10. All 
of these are expressed explicitly in terms of matrix properties. The benefits 
of a matrix (resin) (S m , E m ) relative to a reference matrix (S^q, Em 0 ) are indi- 
cated by a$ ( S ^ / S mo ) for strength and x m (E m /E mo ) for modulus. The mean- 
ing of the parameters C s , C^, C 0 , and C is apparent from the equations in which 
they appear. These parameters are used to include those constituent material 
and ply properties which are independent, or very nearly so, of resin proper- 
ties. Analogous expressions can be obtained for shear and interply delamina- 
tion from the expressions in figure 10 by suitable replacement of variables 
and subscripts. 

From the previous discussion and the expressions in figure 10, the resin 
selection criteria for improvements in individual properties are: 

1. Ac for ply transverse and intralaminar shear strengths, for 
interply delamination and interlaminar shear controlled longitudinal 
compression. 

2. A m for ply transverse modulus, longitudinal tensile and shear 
strengths, and crippl ing-controlled longitudinal compression strength. 
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3. Xc/A m for ply transverse, intralaminar and Interlaminar shear 
energy densities. 

4. *s/*m f° r laminate (composite) strength and "composite tough- 
ness" in general. 

The graphical representation of the "resin selection criteria” is Illus- 
trated in figure 11 where the improvements relative to a reference (state-of- 
the-art), S mo and E^o) matrix for individual and combined properties are 
readily observed. The "resm selection criteria" for improved "composite 
toughness" and structural integrity/durabi 1 1 ty in general can be now simply 
stated as follows: 

"To improve composite toughness, increase the S^/E^ ratio relative 
to S^/Emo with S m increasing at a faster rate than Em". 

It is important to note that S m /E m may be misconstrued as only 
strain to fracture. It is not as is readily observed from figure 12. It is 
also worth noting that the structural integrity dependence on the (S m /E m ) 
ratio is consistent with those obtained by sensitivity analyses in conjunction 

with structural optimization^. 


COMPARISON WITH EXPERIMENTAL DATA AND DISCUSSION 

Available experimental data^» ^0 f 0 r several composite properties made 
from a variety of resins, with stress strain curves as plotted in figure 13 
and in different environments, are compared with the "resin selection criteria" 
described previously. 

Angleplied laminate and unidirectional strengths and moduli are compared 
with A m and a s in figure 14. The correlation for resin-controlled 
properties (transverse and shear) with Ac is excellent. Also the fiber 
controlled strengths (longitudinal flexure), S t nc and S cxx 

correlate with A m . As expected, t^3 fiber-controlled property S^jj does 
not correlate with a™. A very good correlation of resin controlled proper- 
ties (90° (transverse) and shear) with A m only is shown in figure 15. 
Comparisons with elevated temperature data are shown in figure 16. Again the 
correlation of resin controlled properties (transverse and shear) with a s 
is excellent, while there is no correlation of fiber controlled properties 
with A m . 

Various resins, including those from figure 13, are plotted in the resin 
selection criteria space shown in figure 17. It can be seen in this figure 
that 5208 matrix is below the first ply failure boundary and would, therefore, 
be unsuitable for improved fracture toughness. A group of resins, however, 
are identified which will improve significantly the "composite toughness". 

These resins have a s ratios of about 3 and A m of about 2 or A s /A m 
of about 1.5. The delaminated area sustained by a composite under impact 
loading is correlated with the resin selection criteria A s /A m in figure 
18. As can be seen, the correlation is excellent. 

The important conclusion from the correlation results a«. ^ discussion is 
that the "resin selection criteria" derived herein correlates with experimen- 
tal data for a variety of conditions. The resin selection criteria space is a 
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concise graphical means to a priori Identify resins with desirable bulk state 
characteristics which will lead to "tougher" composites. The correlation of 
the criteria with f atigue data Is yet to be determined. However, based on the 
theoretical results, it is anticipated that these resin selection criteria 
should apply to fatigue loadings as well. 


CONCLUSIONS 

Resin selection criteria for tougher composites were der.ved using a for- 
mal methodology consisting of upward integrated and top-down structured 
theories. The criteria account for resin strength and modulus, ply energy 
density, laminate first ply failure, environmental and cyclic load effects. 

The criteria are expressed in resin properties benefits regions where the 
region boundaries are given by simplified equations, or ratios, for resin 
strength and modulus. The resir, selection criteria correlates with experimen- 
tal data for a variety of conditions including laminate strength and stiff- 
ness, elevated temperature effects and resistance to impact. The criteria, 
expressed in a "criteria selection space", provide a formalized direction for 
the a priori selection and the development of resins for "tougher" 
composites. 
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A laminate axial stiffness 

C global damping matrix; laminate stiffness matrix; stress wave speed 

D laminate bending stiffness 

E elastic properties matrix as defined by subscripts; modulus, as de- 

fined by subscripts 

F global force; failure criterion function 

f hygrothermal property degradation factor 

g gravity acceleration 

I Identity Matrix 

i index 

J ply stress influence coefficient as defied by subscripts 

K global stiffness matrix; coupling coefficient in failure criterion 

M global mass matrix; laminate moment as defined by subscripts 

m moisture 

N laminate in-plane force as defined by subscripts 

number of layers in a laminate 
R ply orientation matrix; impacting sphere radius 

5 strength as defined by subscripts 

T temperature 

t thickness as defined by subscripts 

u global displacement 

x,y^ global (structural axes) coordinates 

1,2,3 ply material axes coordinates 

a thermal expansion coefficient as defined by subscripts 

6 moisture expansion coefficient as defined by subscripts 

y global damping matrix proportional ity factor 

a change 

o local indentation 

e strain as defined by subscripts 

e Q global reference plane strain 

5 failure strain as defined by subscripts 

A eigenvalue; resin selection criteria ratio as defined by subscript 

u global curvatures as defined by subscripts 

v Poisson's ratio as defined by subscripts 

p density as defined : subscripts 

12 



0 stress as defined by subscripts 

<p nondimentional function defined by appropriate equations 

u circular frequency 

Subscripts 

C compression 

c composite property 

HTM hygrothermomechanical effect 

1 ply property 

m moisture, hygrothermal effects 

r resin property 

S shear 

s sphere 

T tension, temperature 

xyz respective coordinate directions, properties 

123 ply material axes respective properties 

a T-tension or C-compression 

s T-tension or C-compression 

Matrices 


L 

j 

array, matrix 

{ 

) 

vector, column matrix 

C 

r 1 

matrix inverse 

[ 

j t 

matrix transpose 
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Figure 1. - Defected laminate static fracture data. 
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Figure 2. - Stress concentrations depend significantly on composite moduli. 
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Figure 3. - Governing equations-, micromechanics-elastic continued. 
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Figure 4. - Governing equations: micromechanics-uniaxial strengths. 
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Figure 6. - Simple expression accurately predicts fracture 
stress reduction in perforated plates compared to energy 
method. 
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Figure 7. - Governing equations: micromechanics-hygrothermal properties. 
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Figure 8. - Ply stress influence coefficients. 
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Figure 9. - Ply combined stress failure criterion with hygrothermomechanical effects. 
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Figure 10. - Benefit bounds equations for resin properties. 
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Figure 11. - Resin selection criteria for "tougher" composites. Figure 12. - Resin property regions lor improved 

composite structural integrity. 
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Figure 14. - Resin selection criteria corre- 
lation with measured composite data. 
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Figure 15. - Proposed resin selection criteria 
correlate with measured data. P 0 and E mo 
are reference properties. 


OWOWAL PAGEJS 
OF POOR QUALITY 



3 




<b) Unidirectional composite longitudinal. 


Figuio 16. - Resin selection criteria corre- 
lation with measured composite data at 
elevated temperature. 
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Figure 16. - Concluded. 
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